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Abstract 

We explore in detail the consequences of the CP-violating phases residing in the super- 
symmetric and soft SUSY breaking parameters in the approximation that family flavour 
OhI mixings are ignored. We allow for non-universal boundary conditions and in such a consider- 

ation the model is described by twelve independent CP-violating phases and one angle which 
misaligns the vacuum expectation values (VEVs) of the Higgs scalars. We run two-loop 
renormalization group equations (RGEs), for all parameters involved, including phases, and 
we properly treat the minimization conditions using the one-loop effective potential with CP- 
violating phases included. We show that the two-loop running of phases may induce sizable 
effects for the electric dipole moments (EDMs) that are absent in the one-loop RGE analy- 
\^ i sis. Also important corrections to the EDMs are induced by the Higgs VEVs misalignment 

O^ ' angle which are sizable in the large tan/3 region. Scanning the available parameter space we 

^+ . seek regions compatible with accelerator and cosmological data with emphasis on rapid neu- 

C^ ' tralino annihilations through a Higgs resonance. It is shown that large CP-violating phases, 

OO . as required in Baryogenesis scenarios, can be tuned to obtain agreement with WMAP3 cold 

dark matter constraints, EDMs and all available accelerator data, in extended regions of the 
parameter space which may be accessible to LHC. 



1 Introduction 

Supersymmetry (SUSY) seems to be an indispensable ingredient of Superstring theories and 
super symmetric extensions of the Standard Model (SM) have attracted the interest of physicists 
for more than two decades or so. Supersymmetric models are the only known extensions of the 
SM that are renormalizable field theories, bearing therefore the virtue that radiative corrections 
can be put under control and definite predictions can be made. On the other hand it is well 
known that these models are characterized by many arbitrary parameters, even in their most 
simplified versions, and therefore additional theoretical assumptions have to be invoked to lessen 
their number and build less proliferated models. The minimal supersymmetric standard model 
(MSSM) has the minimal physical content and it needs 124 parameters, a large number indeed. 
These are reduced to much fewer in the supergravity scenarios (mSUGRA), in which universal 
boundary conditions hold, which can be further reduced if additional assumptions are made, as 
for instance absence of generation mixings in the supersymmetric sector at the tree level and/or 
absence of CP-violating phases in the supersymmetric parameters which in the minimal much 
studied versions of SUSY are switched off. For a review on the availability and the description 
of the various models see Ref. [1] and references therein. A thorough account of the parameters 
describing the supersymmetric models and the issue of the CP-violating phases can be found 
in [2]. 

Supersymmetric models with CP-violating phases, other than this occurring in the Cabibbo- 
Kobayashi-Maskawa matrix (CKM), have been extensively studied in the past [3-8]. In the 
mSUGRA models there are only two observable phases, in addition to the CKM phase, which 
are tightly constrained by the EDM data. In more general cases however, when universal 
boundary conditions are abandoned, the number of phases is increased opening new possibilities 
that greatly affect phenomenology. The CP-violating phases residing in the supersymmetric 
parameters produce not only new phenomena, absent in minimal models, but also affect CP- 
conserving quantities, like the mass spectrum for instance, or have large impact on various 
particle processes being therefore of relevance in collider searches [9]. The reconstruction of the 
soft supersymmetric Lagrangian from experimental data, including cases where CP-phases are 
present is addressed in many works [10]. 

The constraints imposed by the EDM data of neutron [11], \dn\ < 6.5 x 10"'^^ e • cm, the 
EDM of electron [12] deducted from measurement of the corresponding EDM of thallium, |(ie| < 
1.6 X 10^^^ e • cm, or diamagnetic atomic systems such as Mercury (Hg-199) [13], \dHg\ < 2.1 x 
10^^^ e-cm, are very tight restricting CP-violating phases to be unnaturaly small. This problem 
is termed as the supersymmetric CP-problem. The phenomenological constraints imposed by 
measurements of the electric dipole moments has been the subject of numerous works [3-6, 14- 
33] u. To forbid overproduction of EDMs one may assume that the masses of superpartners are 
heavy enough beyond the reach of LHC. In other approaches special mechanisms are invoked, 
such as the cancellation mechanism [14-16,18,20], in which contributions of the various Feynman 



^For a thorough review concerning the role of dimension five and six operators affecting the EDMs of atomic 
systems and their link to EDMs of neutron and electron see [32]. For a general review concerning low energy tests 
of the weak interactions including EDMs see [33] . 



graphs involved delicately cancel each other to render EDMs of neutron and electron small within 
their experimental limits. However even in this case the limits imposed by the EDM of Mercury 
atom are hard to satisfy [19,21]. One may cure the situation by lifting the sfermion masses [22] 
but then two-loop contributions are not suppressed [23,24]. The EDMs have been studied 
beyond the one-loop order. The Barr-Zee type [34] two-loop supersymmetric contributions have 
been thoroughly studied [25-29] and yield sizable contributions which do not decouple in the 
limit of heavy sfermion masses. The same holds even for some three- loop contributions arising 
from the running of the renormalization group equations (RGEs) which induce sizable phases 
to the gaugino masses [30]. 

The effect of supersymmetric phases in Higgs sector has been extensively studied [35-41] 
and will not be repeated here; for a review and see [42]. We merely state that their couplings 
to other particles differ from the CP-conserving case and that the CP-odd Higgs mixes with the 
CP-even mass eigenstates. Therefore, except their involvement to EDMs, in principle, they also 
affect neutralino relic densities and thus they are important for cosmological considerations. 

Large phases residing in supersymmetric parameters are welcome for Baryogenesis which 
can occur either through Leptogenesis [43] or through a strong first order electroweak phase 
transition (for reviews see: [44]). Squark and slepton driven Baryogenesis requires a light stop, 
with mass 120 GeV < m^ < rrit, and in conjunction with the fact that the phase transition 
becomes too weak for Higgs masses larger than ~ 120 GeV it leaves a narrow window for 
successful baryon asymmetry [45-48]. Higgsino and Gaugino driven baryogenesis [49,50] is an 
alternative. This effect is resonantly enhanced when the Higgsino mixing parameter /j, is of the 
same order with the gaugino masses Mi^2, [51]. The relevant CP-phases are arg{p,Mi^2) and 
values in the range ~ 10^^ are adequate to produce the observed baryon asymmetry of the 
universe if \fi\ ~ |^i,2|- 

During the last years the WMAP3 [52-54] and SDSS [55] precise determination of the cos- 
mological parameters has stimulated new interest and, in conjunction with the constraints put 
by the accelerator data, it points to a better understanding and more thorough treatment of 
supersymmetric models that violate CP. The Cold Dark Matter (CDM) relic density implied 
by the latest WMAP3 data [54] lies in the range ^cdm^'q = ^-^^^^^oom^ imposing severe con- 
straints on CP-conserving supersymmetric models (for a review see [56] ). The importance of the 
CP-phases in conjunction with Dark Matter (DM) observations has been the subject of many 
works [57]- [58], [50]. In these works the effect of the phases on the neutralino relic abundance 
is discussed observing the constraints put by accelerators in various supersymmetric scenaria. 
For a recent review see [58,59]. 

In this work we focus on mSUGRA-type CP-violating models with minimal flavour violation 
(MFV) and seek regions of the parameter space which are compatible with cosmological and 
EDM constraints and all available accelerator data. In particular we refine the analyses of 
previous works by taking into account: 

i. The two-loop renormalization group running of all phases included which may induce 
sizable effects at low energies having large impact on EDMs. We show that even small 
phases at the unification scale are responsible of inducing large corrections saturating in 



some cases the experimental limits put on EDMs. Such an effect was studied in [30] for 
the particular case of the trilinear soft scalar coupling, whose phase affects the phases of 
the gaugino masses. However other phases with phenomenological interest, notably the 
gluino phase, may influence EDMs, as we shall show, inducing non-vanishing phases for 
the remaining soft parameters. 

ii. Full treatment of the one-loop minimization conditions in the presence of CP-violating 
sources. The one-loop effective potential depends on these phases and correct treatment, 
taking into account all one-loop contributions, shows that a misalignment of the Higgs 
vacuum expectation values (VEVs) is induced which is phenomenologically important. It 
is worth noting that the relative angle between the Higgs VEVs cannot be rotated away, 
even if the Higgs mixing parameter is taken real, by proper U{1) Peccei-Quinn or R- 
symmetries. The appearance of this phase is due to one loop corrections of the effective 
potential. The importance of this misalignment angle for phenomenology has been pointed 
out in previous works [60,61]. Here we refine previous analyses and in conjunction with 
i) we show that it induces phenomena which although in principle small may have large 
impact especially on the EDMs and the relic density of the lightest supesymmetric particle 
(LSP), which is assumed to be the lightest neutralino. 

iii. The effect of strong interaction phases which in principle do not affect electron's EDM at 
one-loop, such as the gluino phase. We show that, due to two-loop RGE running, this 
phase may induce CP-odd invariant phases at low energies, with important consequences 
for EDMs which are absent in the one-loop analysis. Such a phase is phenomenologically 
interesting since it is observable in gluino production and it is known to affect neutralino 
relic densities indirectly through its influence on the bottom mass corrections [62-66]. In 
this work we also argue that the appearance of such a phase may affect the top-down 
approach of the renormalization group running since it may give rise to large corrections 
for the bottom Yukawa couplings having as a consequence the appearance of Landau 
poles [67] in the large tan /3 regime. 

iv. The updated experimental value of the top mass is rrit = 171.4 it 2.1 GeVn, which affects 
the location and shape of the cosmologically allowed funnels which open in regions where 
an LSP pair annihilates through a Higgs resonance [70]. It is known that these funnels are 
very sensitive to the input top (and bottom) quark mass [71]. 

In our treatment we follow a top-down approach, which is the appropriate handling if the 
low energy physics has its origin at Planckian energies. In this sense CP-violating phases are 
not given at the Electroweak scale but are extracted after a two-loop running of the real and 
imaginary parts of all parameters involved, which are inputs at the uniflcation scale. Under these 
circumstances it is interesting to explore where, and under which circumstances, cosmologically 
allowed regions naturally show up, in which the EDM constraints are satisfled, for large phases 
which are relevant for Baryogenesis and other phenomenological issues. In our study we pay 



^ The central value rrit — 171.4 GeV [68] of the top mass has slided down by 0.5 GeV according to more recent 
CDF and DO analyses [69]. Throughout this work the value mt = 171.4 GeV is used. 



special attention to the LSP pair annihilation through a Higgs resonance which is one of the 
leading mechanisms to produce the right amount for neutralino CDM especially in the large 
tan/3 regime. 

In doing so, we have developed a Fortran numerical algorithm, which treats CP-violating 
sources in the MFV scenario in the approximation of neglecting generation mixings from the 
CKM matrix, whose effects are known to be smalo In this scheme soft masses and trilinear 
soft couplings are assumed diagonal in family space. In our study we discuss in detail the issues 
associated with the mass predictions of such models, focusing on a two-loop RGEs running and 
the effect of the CP-phases on quantities that greatly affect the renormalization group flow and 
the extracted low energy parameters. As has been already remarked, in our analysis electroweak 
radiative symmetry breaking is enforced with all one-loop contributions to the effective potential 
duly taken into account when CP-violations are switched on. 

This paper is organized as follows. In Section 2 we give a detailed account of the role of the 
phases in conjunction with the class of models studied in this work. In Section 3 we discuss all 
subtleties associated with the phases that are involved and discuss all pertinent formulas which 
have a large impact on the numerical analysis. In Section 4 we discuss the importance of the 
energy dependence of the phases involved and in Section 5 we discuss the constraints arising 
from EDMs and Cosmology and present our main results. We end up with the conclusions in 
Section 6 and give a summary of our results. 

2 Description of the Model 

In the softly broken super symmetric theory the Lagrangian is split as 

C^ = -^SUSY + -Csoft , 

where all information concerning supersymmetry breaking is encoded in the soft part £soft- 
Proper U(l) and R-transformations of the multiplets involved may be used to eliminate some 
of phases residing in the parameters describing this Lagrangian leaving aside a number which 
cannot be further rotated away. These phases consist an additional set of arbitrary parameters 
with important phenomenological consequences. There are several works describing the situation 
in extensions of the MSSM, in which the presence of such phases is taken into account and 
their phenomenological implications are discussed. In this work we study MFV versions of 
supersymmetric models assuming that soft masses and trilinear couplings residing in £soft are 
family blind. For brevity this class of models we shall coin CPMSSM. In such models the effects 
of the 1st and 2nd generation Yukawa couplings in the running of the renormalization group 
equations (RGE) is small and can be safely neglected from the analysis. 

To start with, it may help to recall the basic field transformations used in order to eliminate 
the redundant phases of the various parameters involved. MSSM based models have C/q(3) x 
C/t/c(3) X C^d(3) and ^7^(3) x JJec'k^) global symmetries acting on the quark and lepton multiplets 



It is customary in literature to term as MFV those models that allow generation mixings only from the CKM 
matrix. 



which can be used to ehminate redundant phases and real parameter from the quark and lepton 
Yukawa coupHngs. At the end six real quark masses, three CKM angles and one CP-violating 
CMK phase are left which are measurable. Thus, in the approximation that the effect of the 
generation mixing is neglected we deal with real Yukawa couplings which are diagonal in family 
space. However in the limit that the Higgs multiplet mixing parameter /x = and the soft 
terms are neglected additional symmetries exist. These are Peccei-Quinn (PQ) global U{l)pQ 
symmetries and R-symmetries Ur{1) under which the bosonic and fermionic components of the 
multiplets involved are transformed by appropriate phase factors. Under U{l)pQ the Higgs 
multiplets Hi^2 have charge 1, all quark and lepton multiplets have charge —1/2 and the vector 
multiplets are neutral, that is carry zero charge. On the other hand the R-charge of these 
multiplets is 1 for the scalar Higgses, 1/2 for the scalar partners of the quarks and leptons 
and zero for the vector bosons. The corresponding Higgsinos, quark and lepton fermions have 
charges by one unit less while gauginos carry R-charge +1. In the literature are used instead the 
symmetries U{l)pQ and U{l)p_pQ. Under the latter the scalar Higgses carry zero charge. A 
useful way to keep track of the changes implemented by these transformations is to assume that 
the Higgs mixing parameter /j,, the Yukawa couplings htfi^r---^ the gaugino masses Ma, the Higgs 
scalar mixing parameter m^, the trilinear couplings At^^r--- and the squark and slepton masses 
squared tti? - are spurion fields transforming in the way shown in Table [U so that the Lagrangian 
C is kept invariant. Phrased in another way, if a U{l)pQ or U{l)p_pQ transformation is carried 
out on the fields involved, i.e. / -^ exp{iQfa)f, the parameters in the transformed Lagrangian 
appear multiplied by exp{—iQa), with Q the charge shown in Table [TJ It is apparent that 
U{l)pQ affects only ^, m^, under which both have the same charge, and C^(l)/j_pn affects fi, the 
gaugino masses and trilinear couplings. None of these affects the Yukawa couplings and the soft 
scalar masses which means that if real in one basis they remain real after U{l)pQ or U{l)p_pQ 
transformations. In a particular basis one exploits the aforementioned transformations to rotate 
away redundant phases of the parameters involved. Which ones is a matter of convention. 
However certain combinations of phases are invariant under these PQ and R-transformations 
and all physical quantities depend on linear combinations of these invariants. In the CPMSSM 
twelve independent invariant combinations can be defined, 

arg{iJ,Mamf) , arg{fj,Aiml*) 

Any other invariant combination is expressed in terms of these. 

As already discussed the Yukawa couplings can be taken real, positive or negative □. Moreover 
if they are real at one scale they remain real at all scales. The reason is that their RGEs are of 
the generic form dh/dlnQ = S h, with S real. The slepton and squark soft masses can be also 
taken real since their imaginary parts are completely decoupled from the theory. In general the 



^ In our convention the VEVs of the Higgses are (-ffi) = vi, {H2) ~ exp{id)v2, with ui,2 real, and trhe 
superpotential, see lA.ll is such that the running masses for the bottom and tau are rrib = —hbV\ and rwr = —h-rV^. 
We assume hb, hr < so that these are positive. Since {H2) is complex the top mass term is htV2 e * ttc + {h.c.) . 
The phase 6 can be absorbed by a chiral rotation resulting to a positive top mass mt — htV2 if ftt > 0. We 
shall adhere to this convention in the following. This chiral rotation modifies other couplings, like giuino-top-stop 
for instance , which will appear to be 6 dependent. 
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Table 1: Parameters and their PQ, R — PQ charges 
Lagrangian part involving the scalar soft masses has the following form, 

^soft = "2 '^Ij ^Uj + (h-c) (1) 

where summation over the squark, slepton and Higgs fields, denoted by Sj, is understood. The 
RGEs of mjj at two- loop order can be found in [72,73]. In the MFV scenario the matrix rriij 
is diagonal and therefore only the real parts of the diagonal elements ma participate in Eq. ([1]) 
and affect physical quantities. Their imaginary parts are decoupled. Thus, only Re{mfj) matter 
whose signs can be either (+) or (— ). In most of the cases the (— ) case leads to potentials 
breaking colour and/or lepton number, or being unstable, and therefore these cases are not 
phenomenologically interesting. 

In order to choose a particular basis to work with, one can exploit the PQ symmetry to make 
m| real. In this basis the invariants under the remaining U{l)j^_pQ symmetry are arg{fiMa) 
and arg{iiAi). Furthermore, by use of U{l)^_pQ one of the phases of /x. Mo, Ai can be rotated 
away and no further rotations are allowed. Which one is rotated away is a matter of choice. For 
instance in mSUGRA, in the presence of CP-violation and with universal boundary conditions, 
one usually chooses to eliminate the phase of common gaugino mass. Mi 12, at the unification 
scale and two phases remain, this of ^ and that of the common trilinear coupling Aq. Since it 
is customary in mSUGRA to work in the basis in which /i, ^0 ^ire complex it is advisable that 
we work in a basis where the phase of /x is not rotated away either, offering the opportunity of 
a direct comparison of CPMSSM with mSUGRA models. Also since in the EDM cancellation 
mechanism, to be described later, we implement rotations of the phases of Mi and M3, in order 
to obtain values for the EDMs of electron and neutron much smaller than their experimental 
bounds, we had better not rotate away these two phases either. 

One should take care of the fact that, in general, phases run with the energy scale. Thus, 
if a phase is set to zero at some energy scale, it may reappear at some other scale due to the 
RGE running. Exception to that are the Yukawa couplings and the /x parameter whose RGEs 
are multiplicative, by real functions, at any loop order. The RGEs of the soft gaugino mass 
parameters are multiplicative at one-loop, but not at the two-loop order, while those of the 
trilinear couplings are not multiplicative already at one-loop order. Therefore, more phases are 
expected to be generated through the RGE running even if at some scale they are vanishing. In 
mSUGRA like models, for instance, three phases are generated for the Mi^2,3 gaugino masses 
at low energy scales, due to the two- loop running, even if some of those are set to zero at the 



unification scale. We shall discuss this issue in detail later on. 

Except the phases associated with the parameters mentioned above, one further phase 
emerges, through loop effects. This is the misalignment angle of the Higgs VEVs which is 
present even if m| is chosen real. To maintain both Higgs VEVs real at one-loop an additional 
rotation of the Higgs fields should be performed, but then this relative phase moves someplaces 
else affecting other parameters. In our approach we take rn^ real at the minimization scale 
Qmin , usually taken to be the average stop mass, by an appropriate PQ rotation. The reality 
of m| simplifies the solution of the minimization conditions as we shall discuss. In addition one 
can exploit the C/"r(l) gauge symmetry of the Lagrangian to redefine fields in such a way that 
the VEV of Hi is real. Thus, one has (Hi) = vi and {H2) = exp{i6) V2, with wi 2 both real. In 
general they can be taken both complex, (Hi) = exp(ixj) Vj, i = 1,2, but only the combination 
Xi + X2 = ^ is observable. In the following we shall work in the basis in which (Hi) is real. In 
this basis it is more appropriate to incorporate 9 with the phase (j)^ of fi and use the combina- 
tions arg{nMa exp(i6')) = (p^ + S,a + 9, and arg{iJ,Ai exp{i9)) = 4>^ + 4>Ai + to express physical 
observables. The reason of doing that is that chargino, neutralino and sfermion masses, as well 
as the EDMs of fermions, the quark chromoelectric moments, and the dimension-6 Weinberg 
operator, depend on the combination (j)^ + 9 rather on (p^ alone [16]. 

In our treatment we take m\ and {Hi) real at the minimization scale, as we discussed earlier, 
but we do not implement a further U{l)j^_pQ transformation to rotate away one of the remaining 
phases at Qmin- The reason is that even if any of these is rotated out at Qmin, it will reappear 
at the unification scale Mqut-, with the exception of the phase of ^, due to the RGE running. 
Therefore we found it more convenient to deal with all thirteen phases, those of fi,Ma,Ai at 
Mqut- The phase of Trig at Mqut and the misalignment angle of the VEVs, 9, are calculable 
and not free parameters. Although legitimate, in this procedure different choices for the input 
phases at Mqut rnay correspond to the same physical situation. Therefore for given SUSY 
inputs we compare the values of arg{iJ,Ma exp{i9)) and arg{fj,Aiex.];){i9)) at Qmin and 9 after 
each run. If they coincide they correspond to the same physical situation and should not be 
double counted. 

Supersymmetric CP-violation has important phenomenological implications and in conjunc- 
tion with other cosmological and experimental constraints deserves further detailed study. For a 
review see [1] and references therein. The presence CP-violating phases affects the cosmological 
predictions for the neutralino relics which are the leading candidates for CDM. They have large 
impact on the bottom mass corrections, Anib, which in turn affects drastically the predicted 
neutralino relic density ^x^o- This effect is more enhanced for large values of tan/?, and for 
small trilinear scalar couplings the important phases are those of the parameter fi and the phase 
^3 of the gaugino mass M3, [62-64]. For large values of the trilinear couplings, at low energy 
scales, their phases play a significant role also. The appearance of phases greatly affects the 
masses of the Higgs bosons as well, and thus they have a large impact on the DM predictions for 
neutralino annihilations near a Higgs resonance region, which is one of the regions favoured by 
the cosmological data . In [66] particular extensions of the mSUGRA, where the magnitudes of 
the soft breaking parameters are universal but their phases are different in general, were found 



to be consistent with EDMs, WMAP3 data, b ^ sj and b — t Yukawa unification in regions of 
the parameter space in which the phases (j)^,^3 have large values. In these considerations the 
cancellation mechanism among the various contributions to EDMs was invoked [14-16,18,20], 
which proved to be a powerful tool to comply with EDM constraints, in the low mo, M^ ^ regime, 
relaxing the stringent constraints imposed on the CP-violating phases although its validity and 
naturalness have been questioned in other works [19]. 

In this work we shall refine the analysis of the CP-violating models, focusing on mSUGRA- 
type models in which all phases are opened up, but magnitudes of the SUSY parameters are 
universal at the unification scale. These models are a subclass of the CPMSSM, and their 
phenomenology has been studied. However the subtleties associated with the two-loop running 
of all phases involved, in conjunction with the delicate treatment of the misalignment angle of 
the Higgs VEVs, which arises from the one-loop effective potential, have not been fully treated. 
The effects arising from such a consideration are important for EDMs and relic densities as we 
shall discuss. 

3 CP-violation in the top-down approach 

In the extended class of supersymmetric models discussed in the previous section, and in order 
to obtain the mass spectrum and phases at low scales one has to run seventy eight (78) RGEs 
for the real and imaginary parts of all quantities involved, including those of the six trilinear 
scalar couplings of the first two generations which are important for the study of the EDMs of 
the light fermions since they are affected by the gaugino masses and the trilinear couplings of 
the third generation species. Due to this dependence non-vanishing phases for these trilinear 
couplings can be developed, even if they are absent at the unification scale, affecting the EDMs 
of the light leptons and quarks with important phenomenological consequences. 

The RGEs for the gauge couplings are certainly real and the Yukawa couplings can be taken 
real as explained in the previous section. The imaginary parts of the soft squark, slepton and 
Higgs masses run with the energy scale but as stated in the previous chapter have no effect on 
physics. Any choice for them leads to the same physical results and for convenience we take them 
zero at the unification scale. The RGEs, up to two-loop order, can be found in [72] and [73] and 
have been adapted to our own notation (see Appendix). These are run from a GUT scale Mqut 
defined to be the point at which the gauge couplings ai^2 unify. We do not enforce unification of 
these with the strong coupling constant although it is left as an option in our numerical code. In 
this procedure we observe radiative electroweak symmetry breaking conditions which are altered 
in the presence of CP-violating sources. The Yukawa couplings of the first two generations have 
little effect and can be neglected from the remaining RGEs in the approximation that the third 
generation Yukawa's dominate. 

In our analysis we follow a top-down approach with input values for the magnitudes and the 
phases of the soft masses and trilinear couplings at the GUT scale. The reasoning behind this 
relies on the fact that these are not known at low energies but they are rather determined from 
the fundamental underlying theory, which describes physics at GUT or Planckian energies. 

As explained in the previous section we shall work in the basis in which by appropriate 



U{l)pQ syminetry the Higgs mixing parameter m^ is real at the minimization scale Qmin which 
we choose to be the average stop masses scale Qi = {nii rrif ) ' . The neutral Higgses develop 
VEVs along the directions {Hi) = vi , (-^2) = '^2 e and it is convenient to shift the neutral 
Higgs components as 

As in the CP-conserving case we write vi = ^ cos /3 , V2 = -^^ sin /? by defining the angle (3. 
The misalignment angle is determined by minimizing the scalar potential and it is vanishing 
at the tree level. The minimization conditions are then given by 

1 , , ,9 „ , m? — TTin tan (3 cos^ 9 

{Ml + Uzz^ - ^ ^ 



2 ' ^ """ ( tan/32 - 1 ) ( c2 + s2 cos2 e ) 

\ cos 8 I mf + 772-2 

2 1 aAy 

rric, smO = —= ;— — . 

^ V2v2 dh 

In the equations above C/j = cos/3, S/3 = sin/3, /SV is the loop corrected scalar potential, and 

_2 2 ^ gAF 

diReHff 

In the first of EqsEJ Mz is the physical (pole) Z-boson mass and Iizz{k) is the transverse Z— 
boson propagator correction at k"^ = M^. Its inclusion is important for a correct numerical 
treatment. Note that the expression on the l.h.s. in the first of Eq. [3] defines the squared of the 
running Z— boson mass, M^ = M^ + Tlzz, which is used to derive the relations of ui,2 with 
the other quantities. When the CP-violating phases are switched off the third of the equations 
above yields a vanishing misalignment angle 9, since its r.h.s. vanishes at the minimization 
point li^Ri = 0. In this case the first two of Eqs. [3] receive their well-known expressions valid 
in the CP-conserving case. In our treatment all one-loop corrections to the effective potential 
/SV have been calculated allowing for CP-violating sources in the field dependent masses of all 
SUSY sectors involved, i.e. sfermions, charginos, neutralinos and Higgses. 

In our numerical approach we have tan /3 as input and the value of 77i| at the minimization 
scale is determined by the second of Eqs. O The magnitude of ^ parameter is output determined 
by the first of Eqs. [3l Recall that m^ = |/x| -|- rn|^., z = 1, 2, with mui denoting the soft Higgs 
masses. In the presence of CP-violating sources the phase 6 is non-vanishing, because of loop 
corrections to the effective potential, and it is determined from the third of Eqs. [3]( see [60,61]). 
Thus, it is expected to be small but its impact on the electric dipole moments of known species 
may be sizable. Since this is a one-loop quantity one can consistently solve Eq. [3] by putting 
6* = within /SV . 

The fields -Ri, 2 1-^1,2 in Eq. [2] are not mass eigenstates. A linear combination of /i^2i namely 
I2 = —cpli + SJ3I2 is the Goldstone mode and the orthogonal to it, I[ = spli + 0^/2, gets mixed 



with Ri2 through a 3 x 3 mass matrix. When CP violating effects are absent this matrix does 
not allow mixing of I[ with Ri^2- In this case I[ is the pseudoscalar, CP-odd, mass eigenstate. 
The other modes Ri^2 do get mixed and they must be rotated to yield the heavy and light 
CP-even mass eigenstates . 

The corrections to the masses of the third generation, through which the corresponding 
Yukawa couplings are read, are very important and affect the numerical treatment. The presence 
of CP violating sources affects these corrections substantially. The supersymmetric corrections 
to the bottom mass are sizeable for large values of tan /3 and should be duly taken into account 
in the analysis. These give rise to large corrections to the bottom Yukawa coupling [74-76] given 
by 

\ rib\ - (1 + ^sUSY) ■ (4) 

Vi 



Throughout with hat we denote quantities in the DR scheme. In the equation above the SUSY 

corrections are denoted by ^'susY ^^^ they are resumed according to the scheme presented 

in [75]. In this equation rhf {Mz) is the Standard Model DR value of the bottom running 

mass at the scale Mz- Its value is calculated by running the S'C/c(3) x C/em(l) RGEsu for masses 

and couplings in the MS scheme, from the bottom mass m^ (rrib) = 4.25 it 0.15 GeV [78], as 

determined in lattice calculations, up to the scale Mz- Its MS value at Mz is subsequently 

converted to its DR value using well-known expressions. Therefore from Eq. (j3|) the value of hf, 

can be extracted which is needed to run the RGBs from Mz to the GUT scale. The corrections 

involved within A^^^y are very important and are discussed below. 

The leading supersymmetric QCD, sbottom-gluino, and Electroweak (EW), stop-charginos, 

contributions to A^j^^y are given by 

2a ~ ~ 

^SUSY = -ir^Mg{\fi\ta.npcos{^3 + ^^ + 0) + \fi\\Ab\cos{^3-A))G{bi,b2,Mg) 
ovr 

- TTni 1^1 ( lMl + l^t|tan/3 cos(</.^ + 0* + 0)) G(ti, ts, |Ai|) ■ (5) 

|Ab^t|,|;u| are the magnitudes of Ah^t,fJ' and (/)b^t,0^ their phases. 9 is the misalignment angle 

between the Higgs VEVs and tilded quantities denote sbottom, stop masses |j. In this expression 

we have neglected electroweak (EW) mixings of the stops, and also sbottoms, and the mass of 

the chargino is put to fi. More refined expression which take into account the mixings can be 

found in ref [76]. An analogous treatment holds for the corrections to the tau mass as well. 

However their effect is small due to the absence of QCD corrections at the one-loop level. When 

the magnitudes of the trilinear couplings involved in Eq. [5] are small and the angle 6, which is 

anyway small, are neglected then this equation receives a much simpler form 

2a ~ ~ 

^SUSY = ^M~g |/i| tan/3 cos(6 + </>/.) G{bi,b2,Mg) 



3tt 



h? 



167r2 



* Vl'G'(ii,t2,|/u|). (6) 



^ In our treatment we use 3-loop QCD and two-loop QED RGEs for the strong and electric couplings and 
two- loop for the running masses of the bottom and tau. Four- loop, ©(of), QCD contributions to beta functions 
and quark anomalous dimensions are available [77] but affect little our analysis. 

^The function G{a,b,c) is identical to I{a,b,c) defined in Eq. 7 in [75]. 
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The functions G appearing in Eq. [6] are positive in a large portion of the parametric space and 
therefore if cos (^3 + (f)^) is negative the corrections to ^\usy ™^y turn out to be negative and 
sizable, in the large tan (3 regime. In this case the bottom Yukawa coupling of Eq. H] gets large 
and a Landau pole may develop. Therefore in the top-down approach, depending on the inputs, 
the approach to the large tan (3 regime is not guaranteed. This observation is important for the 
mechanism of neutralino annihilation through a Higgs resonance, which opens for large values 
of tan (3. The importance of this will be discussed later. 

The Yukawa coupling of the top quark is large and it is determined from the experimentally 
measured top pole mass in the following way. For the top quark the relation between its pole 
and running mass, including the dominant QCD and the supersymmetric gluino-stop corrections 
is given by 

Mf' = m,{Mf''){l + ^'susY)l{l-^QCD). (7) 



The pole mass is scheme independent so that the r.h.s can be calculated in either MS or DR 
scheme. We prefer to employ the first scheme since the QCD corrections ^qcd have a simple 
form which at two-loops is given by 

Aqcd = ^a« + l.llas^ . (8) 

ovr 

In Eq. [8] the strong coupling constant is meant at Mf° ^ and its running from any lower scale 
to the pole top mass is done using the five quark flavour RGEs. Note that in Eq. [7] the 
QCD corrections have been resummed. The MS strong coupling ag appearing in the above 
expressions is different from the corresponding DR strong coupling usually denoted by 03. The 
relation between these two will be given in the sequel. 

The gluino-stop corrections appearing in Eq. [7] are given by 

^*SUSY = ^[-B,{0,Mg,mi^)-B,iO,Mg,mi^) 



Sir 



^ sin(20t) cos^ { Bo{0, M~g,m i^) - Bo{0, Mg,m i^) ) 



mt{Mr'') 

In this 6t is the angle diagonalizing the stop mass matrix, ^ = ^3 + 4> + and rn i are the 
stop masses |I|. This generalizes the results of ref. [79] in the case that the supersymmetric 
parameters are complex. The functions Bqi are defined as in [79]. A minus sign difference with 
the results of that reference, occurring in the case of the absence of CP violations, ^ = 0, is due 
to the slightly difference notation used here. Note that these corrections have the same form in 
both MS and DR schemes since in the gluino-stop loop scalar particles are exchanged and no 
traces of gamma matrices are involved. Their difference in the two schemes is therefore small, 
two- loop, due to the fact that the couplings and running masses appearing already differ at 



^ The diagonalizing matrix K is defined by KM^K^ = diag {m?i , m?i ) with the matrix elements of K given 
by Kii = K22 = cos 9t, K12 = —K21 = e^ sindt. The definition of the matrix K we adopt is consistent with 
K tending to the unit matrix if the off-diagonal elements of the squark mass matrix are switched off. Thus, the 
subscripts in m, j do not specify the ordering of their heaviness. 
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one-loop in the two schemes. More refined relations including the subdominant EW corrections 
are given in [76]. The EW supersymmetric corrections are however small and they correct the 
approximate result by less than 1% which becomes even less when the stop masses are of the 
order of 1 TeV. Therefore the MS value for the top Yukawa coupling is given by 

(9) 



t y^^^t ) - — ;:— 1 _L At 

SUSY 



V2 1 + A 



MS 



and its DR value needed to run the corresponding RGEs is provided by the usual conversion 
formula [80] 

hfS^^pole. ^ hf^mr^') ( 1 + ^ + ^^ + ^ ) . (10) 

* ^ * ^ * ^ * ^ ^ Svr 167r SOvr ^ ^ ^ 

Besides the input Yukawa couplings, we need the values of the gauge couplings at the unifi- 
cation scale Mgut defined to be at the point where the gauge couplings ai_2 meet. These are 
determined by their values at the scale Mz in terms of the fine structure constant ao) the value 
of the Fermi coupling constant Gp and the physical Z-boson mass Mz- These are related to 
"^1,2 (^z) in the way prescribed in [79]. These are run up by two- loop RGEs and the unification 
scale Mgut is determined at the point where these intercept. Naive gauge coupling unification 
would entail to putting the strong coupling a^ equal to ai^2 at Mgut- Instead of doing this we 
follow the alternative, followed in [79], according to which the MS value of the strong coupling, 
denoted by ag, coincides with the one experimentally measured. Its relation to 03 at the scale 
Mz is given by 

as{Mz) = as{Mz)/{l - Aas) 

where Aas includes supersymmetric threshold corrections and constants associated with passing 
from MS to DR scheme. This determines 03. The Yukawa and gauge couplings determined at 
lower scales in the way prescribed earlier are run up at Mgut to determine their values at the 
unification scale. This is done iteratively until a certain convergence is achieved. In our code 
the Higgsino and Higgs mixing parameters fi, m^ are outputs, and the convergence criteria in 
our numerical code are tailored to monitor these parameters in each iteration, which in most 
of the cases have the slowest convergence from the other parameters involved. In determining 
these we use the full one loop effective potential with the leading two-loop QCD correction taken 
into account. The minimization conditions are considered at an average stop scale as described 
earlier and the Higgs mixing parameter m\ is taken real without loss of generality as we have 
discussed. At any other scale, including the unification, this is certainly complex and its phase 
is extracted from the running of the RGEs from the minimization scale. 

Concerning the mass spectrum of SUSY particles the gluino physical mass is overwhelmed by 
large QCD corrections which affect the numerical analysis and should be discussed. These are 
due to both SM gluon exchanges and supersymmetric corrections due to exchanges of fermions 
and their corresponding squarks. The relation between the physical, M^, and the soft gluino 
mass, M3, is found to be 

M- = \M^^ (11) 

' (l-3a3(Q)(5 + 61n(Q/|M3(Q)|)-5(Q))/(47r) ^ ^ 
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where 03 is the DR value of the strong couphng and S{Q) the squark contribution given by 

1 3 

S{Q) = - Y,{2q, + n, + di) . (12) 

"^ i=i 

In S{Q) Qi denotes the contribution of the squark doublets, accommodating the left-handed up 
and down squarks, and Ui,di are those of the corresponding right-handed squarks. In both cases 
the index i runs over the colour. For qi, and the same holds for Ui,di, 

I M^ 1 (r - if 1 

where Mj = max{rhi, \M^\ ), r = IM3I /rfif. rhi is the squark soft mass in each case. We have 
neglected EW mixing effects which have little effect on this formula. This generalizes the result 
of [79], in which a common squark mass is used, and it is a handy expression to use avoiding 
the complexities of other calculations. More refined two-loop corrections in terms of several 
two-point functions are presented in [81,82]. In that reference it is shown that the two-loop 
QCD corrections are small. 

With the above we end the discussion concerning the treatment of the RGEs by giving an 
outline of the salient features which have a large impact on our numerical analysis in the CP- 
violating case. The boundary conditions employed for the couplings at the unification scale 
were discussed earlier. In our code the soft masses and trilinear scalar masses we treat in the 
most general case allowing for non-universal boundary conditions for their magnitudes and their 
phases, as already discussed, in the approximation of neglecting flavour violating interactions in 
the Lagrangian. Although the setup is to cover the most general case in this respect we shall 
only discuss particular cases which are of phenomenological and theoretical interest, and focus 
our attention mainly on models with universal boundary conditions for the magnitudes of the 
SUSY breaking parameters involved. 

4 Running of the CP-violating phases with the energy scale 

The case of non-vanishing phases has features that need be further discussed when the two- 
loop RGEs are used for their evolution. Their values are inputs which can be set at either the 
unification scale, Mgut, or the average stop scale Qf at which we minimize the one-loop effective 
potential. In our treatment we work in the basis in which the phase of mg is put to zero at 
Qj and the remaining phases we take as inputs at Mgut, as we have already stated. Evidently 
the phases of all parameters at any other scale are determined by the RGE running. The phase 
of the /U is very important since it affects the mass spectrum of the charginos, neutralinos and 
sfermion and it explicitly appears within their corresponding mass matrices. Besides, its phase 
has a large impact on the radiative corrections of the bottom and top Yukawa couplings and 
therefore the numerical procedure is very sensitive to its input value. Last, but not least, the 
phase of fi affects the EDMs of Hg, neutron and electron. A large phase for n may be in accord 
with Baryogenesis but its value is severely constrained by the experimental EDM bound of the 
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electron in mSUGRA models in which only the common gaugino phase and the phase of /i are 
present. 

The RGEs of the parameter /i, the gaugino masses and the electron's trilinear coupling, 
which also affect the light fermion dipole moments, are as follows 

„ „ = -— (-3a2 - ^ai + 3/i? + 3/i? + hi + two - loop) 
alnQ 47r 5 

dMi 

— — - = PiMi + two -loop, i = l,2,3 
alnQ 

" - (-602^2 - — aiMi + 6hlAb + 2hlAr) + two - loop . (14) 



dlnQ Att 5 

For lack of space we do not present the RGEs of the remaining trilinear couplings of the 1st 
and 2nd generations and we also avoid presenting explicitly the two- loop contributions. The 
Pi in the RGEs for the gaugino masses are the one loop beta function for the gauge couplings 
aj. Since we work in the MFV scenario we have neclected generation mixing terms. In this 
approximation the trilinear couplings of A^, as well as the remaining trilinear couplings of the 
first two generations, do not have any influence on the remaining parameters although they 
depend on them. In practice this means that one can first solve the RGEs for the rest of the 
parameters and subsequently determine the trilinear couplings of the first two generations. 

The RGEs for the fi parameter has the form dfi/dlnQ = ^ S where S is real. Due to this the 
phase of // does not get renormalized with the scale. Its value at any scale is the same with its 
value at Mqut independently of the other phases. The phases of the soft gaugino masses exhibit 
a slightly different behaviour. Their one-loop RGEs are multiplicative as in the fj, case so that 
they do not get renormalized at this loop order independently of the remaining phases. However 
this does not hold at the two-loop order. This essentially means that one should expect small 
renormalization of their phases as we run down from the unification to lower scales and vice 
versa. Thus, even if their phases at the unification scale are set to zero, small phases are induced 
at low energies, if the phase of /x or the trilinear couplings are non- vanishing at Mqut- Although 
small this phenomenon may have dramatic consequences for the EDMs since even small phases 
may produce large EDMs which are constrained by the data. This effect is more enhanced for 
the phases of the trilinear couplings which are considerably renormalized since their RGEs are 
not multiplicative, already at the one-loop, unlike the parameter // and the gaugino masses. For 
instance in the case of A^ its one-loop RGE shows a dependence on the gaugino masses and 
the trilinear couplings of bottom and tau, as is obvious from the third of Eqs. O Thus, a 
non- vanishing phase, for at least one of them, at the unification scale may yield a non-vanishing 
phase at the EW scale for Af, which in turn affects the EDM of the electron. The same holds 
for the remaining trilinear couplings. This digression shows how important is to consider the 
running of the phases as we do in our approach. 

As a preview of the impact of the two-loop RGE corrections to the phases, and consequently 
on EDMs, consider the particularly interesting case arising when the gluino phase ^3 is the 
only non- vanishing phase at Mqjjt- As is obvious from the above RGEs the phases ^1,^2 of 
Ml , M2 as well as this of A^ are not affected by ^3 at one-loop running. Therefore at this 
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Ml/2 = 900 , mo = 800 , ^o = 500 , m* = 171.4 , rub = 4.25 



6 = 27r/10 



tan (3 = 50 
taii/3 = 10 



arg{Mi) arg{M2) arg{Ae) 6 10^'' x d, 






-0.0144 -0.0198 -0.1735 0.0465 +1.0845 
-0.0149 -0.0209 -0.0079 0.0006 -0.1654 



Table 2: The induced phases, at Mz, of Mi^2, ^e, the misahgnment angle 9 and the value for the 
electron EDM when only the gluino phase ^3 = tt/5 is switched on at the unification scale, for 
tan/3 = 50, 10 respectively. The inputs (in GeV) are shown on the top of Table. Mi/2,mQ,AQ 
are the magnitudes of the common gaugino and scalar masses and trilinear couplings at Mqut- 

loop order the electron EDM bound does not depend on ^3. However at the two- loop order 
the phases for Mi,M2,Ae do depend on ^3 affecting the EDM, de, of the electron. In some 
cases this dependence may have important consequences inducing corrections to EDMs that 
are comparable to those induced by the Barr-Zee type two-loop contributions which are known 
to be sizable for large A and tan/] > 30. If, for instance, .^3 is the only non-vanishing phase 
switched on at the unification scale the induced phases for Mi , M2 , A^, may result to values for 
electron's EDM that can even saturate the experimental limits put on de- A typical case is 
shown in Table [21 for tan /3 = 50 and tan /? = 10 respectively. The remaining inputs are as 
displayed in Table. Throughout this paper, if not otherwise stated, Mi/2,mo and Aq will denote 
the magnitudes of the common soft gaugino and scalar masses and the common trilinear scalar 
coupling respectively. For ^3 = 27r/10 at Mqut the phases of Mi^2 at the EW scale are ~ 10~^. 
The phase of A^ is 4>Ae ~ 10~^(10"^) for tan/3 = 50(10). For completeness we also display the 
value of the calculated misalignment angle which also affects EDMs. It is seen that values for 
de which are larger than the experimental limits quoted in literature, \de\ < 1.6 x 10"^''', are 
induced for these particular inputs. However this is a generic feature valid in a large regions of 
the tuq — Mii2 plane. In Fig. [1] we display the ratio \de/dT^\ of the predicted electron's edm 
to its experimental bound as function of the gluino angle ^3, which is assumed to be the only 
non- vanishing input CP-violating input at Mqjjt, for two different values of tan/3 = 20, 40. 
The remaining inputs are displayed in the figure. This ratio should be less than unity for the 
experimental bound to be observed. A strong dependence on the angle ^3 is seen which is absent 
at the one- loop RGE running. This restricts the allowed .^3 values to be in the vicinity of 0, ibvr 
for tan/3 = 20. Increasing to tan/3 = 40 the dependence of de is still strong, but slightly milder, 
and the allowed range for ^3 gets broadened. In this case one observes that large tan/3 values 
allow for non-trivial gluino phases which are welcome since they can be used in order to lower 
the neutron's edm by using the cancellation mechanism as we shall discuss. 

The dependence of the Mi, M2, Ae phases on ^3 may also affect the cancellation mechanism 
if phases are inputs at the GUT scale. Rotating the phase ^1 till de becomes vanishingly small, 
to comply with its experimental bounds, a subsequent rotation of .^3 at Mgut, in order to 
make neutron's EDM small within its experimental limits as prescribed in [20], will shift the 
initially found ^1 invalidating the cancellation between neutralino and chargino contributions in 
de- Certainly this is not the case when this procedure is implemented with phases given at the 
EW scale. Therefore the determination of cancelling phases in the top-bottom approach poses 
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(a) 
i,„„=800, 600 ,.. tan|S=20 
11,^=171.4,4.25' 
\= 100 



(b) 

m„j^=800, 500 
A„= 100 
t|, = 0.0 
m, = 171.4,4.25 



tan|3=40 



Vi 



Un 



Figure 1: The ratio {de/df^^l of the predicted electron's edm to its experimental bound as 
function of the gluino phase ^3 , for two different values of tan /? = 20 (left panel) and tan /? = 40 
(right panel). The remaining inputs are displayed in the figure. 

difficulties not encountered in the one-loop RGE running. In addition, even if proper phases 
are found at Mqut by the cancellation mechanism, so that electron and neutron EDM become 
small for some particular SUSY inputs, it is difficult to delineate regions by merely rescaling 
the SUSY parameters as prescribed in [20]. The reason is that due to their RGE running the 
induced phases at low energies depend on SUSY inputs and a rescaling dislocates the values 
of the low energy phases by little amounts but enough to make the cancellation invalid. We 
are therefore arguing that the cancellation mechanism is best suited for a bottom-up approach 
although it may still be a powerful tool to locate regions compatible with EDMs and all other 
data when CP-violating phases are switched on as we shall discuss. 

The role of the misalignment angle 9 to the EDMs should not be passed unnoticed. This 
is measurable and cannot be rotated away [1,60,61]. Its value enters and affects various phys- 
ical quantities. The neutralino, chargino, squark and slepton masses depend on this through 
the combination arg{ij,) + 9 as has been already stated. It also affects the Higgs decays to bb 
enhancing the widths for these decays [60] which has important implications for the cosmolog- 
ically acceptable regions in which LSP pair annihilation takes place near a Higgs resonance. 
This mechanism depends sensitively on the corresponding widths. It has also impact on EDMs 
especially in the large tan/3 regime. In Fig. [21 for some particular inputs, we display the values 
of the angle 9 and the electron dipole moment, de, with and without the inclusion of 9 in its 
calculation. The angle 9 takes values from ~ 10~^ to ~ 10~^ for tan/3 in the region 5 — 50. On 
the left panel the gluino phase has been taken vanishing at the GUT scale and the difference in 
dgreaches 50% for tan f3 = 50. On the right panel in addition the gluino phase is switched on and 
the difference gets much larger. Therefore the misalignment angle produces large effects when 
tan (3 is large, which are further augmented if in addition the gluino phase is large at the unifi- 
cation scale. Since 6 has large impact on EDMs, in particular regions of the parameter space, 
it influences the cancellation mechanism as well, in the large tan /3 region, especially when large 
values for the gluino phase are required to make neutron's EDM small within its experimental 
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Figure 2: The misalignment angle 9 (solid line) and the electron dipole moments, with 9 included 
(dashed-double dotted line) and with 9 set to zero (dotted line), for the inputs shown on the 
figure. Masses are given in GeV and Mi 2,3,^710,^0 refer to the magnitudes of the corresponding 
parameters. On the left (right) panel the M3 phase is set to (vr/S). 
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Figure 3: One of the heavy neutral Higgs mass (solid line) and the value of twice the LSP 
neutralino mass (dash-dotted line) as functions of tan/3 for the inputs shown in the figure. 
Above tan/3 ~ 41.0 a Landau pole is developed. 

bounds. 

As already stated in this section, the phases of the Higgsino parameter // and that of the soft 
gluino mass ^3 affect the analysis a great deal. The first, if large, affects the EDM of electron 
which imposes the stringent constraint on the phase (/>^ of /i, while both affect the corrections 
to the bottom mass especially for large values of tanfi having a large impact on the DM relic 
density. This is clearly seen in Eq. [5] or its more simplified form [6] where the corrections to 
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bottom mass are encoded in. From Eq. [6]we see that if (j)^, ^3 are such that cos{(j)^ + ^3) < 0, at 
low energies, these corrections, depending on inputs, can be large and negative, in the large tan/? 
regime. Therefore in view of Eq. [Hthey may yield large values for the bottom Yukawa coupling. 
In this case one should be prepared to encounter the appearance of Landau poles and the top- 
down approach can not be handled perturbatively. Therefore the link between low energy and 
GUT scale physics is questioned in this case. This behaviour imposes a severe obstacle when 
large phases are sought in conjunction with large tan/3 values which is the requirement for the 
LSP annihilation through a Higgs resonance. In Fig. [3] we present such a situation for the inputs 
displayed on the figure. In this figure we observe the evolution of one of the heavy neutral 
Higgses mass, rriH^, and the double of the LSP neutralino mass, 2 misp, as functions of tan/3. 
The Higgs mass tends to 2 misp as tan /? increases and would eventually catch the 2 misP line 
signaling approach to a point where LSP pair annihilation through a Higgs resonance dominates 
the relic density. However it is shown that this is abruptly stopped due to the development of a 
Landau pole at around tan/3 ~ 41. This effect, not considered in previous analyses, may exclude 
particular points in the parameter space and shrink the allowed funnel regions. 

5 Cosmologically and EDMs allowed domains 

In the previous section we gave an account of the salient features of MFV models when CP is 
violated. Our principal aim in this section is to explore regions of the parameter space in which 
the LSP neutralino relic density is within the stringent limits put by WMAP3, satisfying at the 
same time the EDMs and all available accelerator constraints. We focus our analysis on regions 
where the LSP neutralinos are paired annihilated through a Higgs resonance which is one of the 
prominent mechanisms in CP-conserving models, in order to obtain acceptable relic densities, 
although other regions of interest will be also explored. Regions of the parameter space in which 
such a process is feasible, when CP is conserved, have the shape of funnels lying on each side of 
the line on which 2mLSp/^A is unity and this occurs for large values of tan /3. In the presence of 
CP violation, and especially for large values of the phases which we are interested in, the shape 
and the location of the funnels are different due to the impact of the phases on the corrections 
to bottom mass as has been discussed in the previous section. There are two cases which one 
should be interested in: 

A. Cases where EDMs are naturally suppressed, without invoking any special mechanism, to 
lower the values of EDMs to acceptable levels having at the same time some of the phases 
large. These regions correspond to large mQ,Mii2 > ©(few TeV). Certainly the EDMs 
can be suppressed if all phases are very tiny but this case is not physically interesting. 

B. Cases where EDMs are suppressed due to cancellations among the various contributions. 
This requires a tuning of the phases at low energies that have a large impact on EDMs. 

In general EDM constraints exclude a large portion in the tjiq, Mi 12 plane allowing large values of 
rriQ or M1/2 so that the EDMs are naturally suppressed. Thus, depending on inputs they cut off 
substantial part, or all, of the cosmologically allowed neutralino and stau coannihilation [83] tale 
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Figure 4: The cosmologically allowed region (shaded light green) for mt = 169.3 GeV (left 
panel) and mt = 173.5 GeV (right panel), when all the phases are switched off, for tan/3 = 50. 
The magnitude of the common trilinear scalar coupling is taken Aq = 100 GeV. The remaining 
inputs are shown on the figure. The solid black line on the left of the figure is the chargino mass 
bound rric > 105 GeV. The red dashed (blue dashed-dotted) line indicates the Higgs mass bound 
114 GeV (115 GeV). The hatched area at the left-top designates the no-electroweak symmetry 
breaking region. At the bottom the shaded region is excluded since there the stau is the LSP. 

and focus point region [84] as well. At the same time they cut part of all of the cosmologically 
allowed region in which neutralinos annihilate through a Higgs resonance. These regions have 
the shape of funnels, whose location and form is sensitive to the top and bottom masses and 
may occupy regions allowed by the EDM constraints, if they happen to span large mo,Mi/2 
values. The sensitivity with the top mass is shown in Fig. HI On the left panel and for the 
inputs shown we display the cosmologically allowed region for m-i = 169.3 GeV, the lowest 
allowed by the experimental data for the top mass. On the right panel the same figure is 
shown with rrit = 173.5 GeV which is the upper experimental limit. In both cases all phases 
are switched off but this behaviour holds in CP-violating cases as well. In the first case the 
cosmologically allowed regions, which follow the MHiggs/^misp = 1 curve, are not so peaked 
and occupy regions characterized by M^/2 < 900 GeV. Increasing mt, on the right panel, the 
location of the line MHiggs/'^T^LSP^ which controls the location and the shape of the rapid Higgs 
annihilation region, turns to the right dragging with it the cosmologically allowed region to 
higher M^ /2 becoming more pronounced and extended. This behaviour is due to the sensitivity 
of the Higgs mass spectrum with mt- The larger the top mass is the sharper the shape of the 
funnel, which extends towards large mo,M^/2 values, and larger the possibility of overlapping 
with EDM allowed regions. For the bottom mass the tendency is rather opposite and it is low 
values of m^ that favour the formation of sharp cosmologically allowed funnels. 

Prom the previous discussion it becomes evident that regions allowed by both EDM and DM 
constraints are easier to find for the highest allowed values of the top mass and values of the 
phases that make the running bottom mass minimum. In previous works [66] the values of the 
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Figure 5: The effect of the gluino phase alone for tan/3 = 45. On the right (left) panel we take 
^3 = 0.1 TT (^3 = 0). The remaining phases are zero and Aq = 100 GeV. On the right panel 
all displayed region is excluded by neutron and Hg EDM bounds. The allowed region by the 
electron EDM bound, de = 1.6 x 10~^^e • cm, lies between the blue short-dashed lines. The rest 
of the curves and shaded regions are as in Fig. [H 



top quark was taken as large as 178 GeV in agreement with the experimental values quoted 
at that time. In view of the new experimental values of mt, according to which mt is lowered 
by about ~ 5 GeV, and due to the sensitivity on the top mass of the cosmologically funnels 
for LSP annihilation through a Higgs resonance , this picture may be distorted. It should be 
also remarked that in our analysis the situation is different from that encountered in [66] where 
Yukawa unification is enforced. In that case Yukawa unification entails to a different bottom 
mass value at each point of the mo — M^ji plane and therefore due to the sensitivity on m;, our 
findings cannot be directly compared to those of [66]. 

Before embarking on presenting our results we remark that for the calculation of the electric 
dipole moments and the Higgs masses we use the FeynHiggs-2.5.1 code, [85]. However since (^2) 
is not real, in order to implement the effect of the misalignment angle 6, one has to replace the 
phase (j)^ of // by (/>^ + 9. For the electric dipole moments of the known species the correctness 
of this we have also checked numerically by comparing the outputs of our numerical routines 
against to those returned by FeynHiggs. For the Higgs masses, in all cases studied, the Higgs mass 
spectrum obtained was close to that obtained by using the effective potential, with an accuracy 
2 — 5%. The latter is known to be less accurate, since the wave function renormalization effects 
are not counted for. Therefore in this work we use the Higgs masses as returned by FeynHiggs 
and this comparison serves as a further check of the correctness of our treatment. Concerning 
the experimental limits put on Higgs masses, the ratio ^^ = {ghzz/gh¥z) °^ ^^^ light Higgs 
boson h coupling to Z, ghzz, to that of the corresponding SM coupling 5^^^ is considered for 
each case studied. This was found to be very close to unity, that is the light Higgs acts like a 
SM Higgs boson, and therefore the LEP2 limit m/i > 114.5 GeV applies. This is expected in 
the constrained model studied in this work since we are within the decoupling region. 



20 



In case A in order to locate regions compatible with all available data the funnel must 
be extended towards high mo,M^/2 entering regions in which EDMs are naturally suppressed. 
However, we have found that this cannot occur since EDM bounds require very high mQ,Mi/2 
not overlapping with the cosmologically allowed funnel regions, unless, the CP-violating phases 
are very small. In order to conceive the picture, and for demonstrating the importance of a 
non-zero gluino phase, on the left panel of Fig. [5] we display the cosmologically allowed regions, 
having the shape of funnels, when CP-violating phases are absent. The magnitude of the common 
trilinear scalar coupling is taken Aq = 100 GeV. The remaining inputs are shown on the figure. 
The hatched area on the left is not allowed due to the absence of EW symmetry breaking 
there. The solid black line on the left of the figure is the chargino mass bound mg > 105 GeV. 
The dashed line (in red) sets the Higgs mass boundary line niHiggs > 114 GeV, while to the 
right of it the dashed-dotted line (in blue) designates the bound on the Higgs mass reported 
by DO [69], niHiggs > 115 GeV. On the right panel for the same inputs we switch on the 
gluino phase ^3 = 0.1 vr. The funnel is slightly deformed with its top end approaching the point 
mo = Ml/2 = 2000 GeV. However all displayed region is excluded by neutron and Hg EDM 
bounds and it is not of physical relevance. The electron EDM, dg, is also affected and allows the 
region confined between the short-dashed lines (in blue). This figure demonstrates in a clear way 
the impact of the gluino phase on de by the two-loop RGE running of the phases, as discussed 
in the previous section. 

The other option opened, case B, is to implement the cancellation mechanism according 
to which phases are chosen so that the various contributions to EDM cancel each other. This 
option does not require high values for mo,Mi/2. Starting for such a point in the parameter 
space one can draw trajectories, by merely rescaling the values of rriQ and M1/2, in which EDMs 
are in agreement with the experimental bounds put on them [20]. These trajectories may 
eventually overlap with the cosmologically allowed portions and thus delineate regions allowed 
by both EDM and cosmological constraints. In the previous section we discussed that this is 
rather difficult to be accomplished in the top-down approach due to the two-loop running of 
the phases and the interplay among these. Certainly the phases can be tuned at the EW scale 
to cancel separate contributions in electron and neutron dipole moments. If such values are 
obtained at the EW scale, when one is subsequently trying to find extended regions by rescaling 
as niQ, Ml/2 ~^ -^ "^Oi -^ M1/2, the RGE evolution of these phases at the GUT scale yields values 
that depend on the parameter A. Therefore their values at Mqut are fine tuned since they are 
different for different pairs of mQ,Mi/2 that are related by a rescaling factor. Therefore in the 
top-down approach it is difficult to find extended regions in the mo, M1/2 plane by just rescaling 
the supersymmetry breaking parameters as prescribed before. 

In our approach, for given inputs for the supersymmetry breaking parameters, we prefer 
to implement this mechanism by rotating the input phases Ci,2,3 at the GUT scale until an 
acceptable point is found respecting the EDM limits on electron and neutron. Subsequently 
keeping fixed the values of the phases we vary m,o,Mi/2 to delineate regions compatible with 
all available data. Alternatively one can vary the phases around the values for which small 
EDMs for the electron and neutron were obtained, keeping the remaining inputs fixed, in order 
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Figure 6: The EDMs and relic density contours for rriQ = M 
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tan/? = 50, and ^2 = 0.007 vr on the (^3,^1) plane. The remaining phases are zero at the GUT 
scale. The blue dashed (orange long dashed-dotted) curves designates the de (dn) acceptable 
region. The allowed region by dng is this between the magenta dotted lines, which does not 
overlap with the grey region allowed by electron and neutron EDMs. The contours of constant 
neutralino relic density are shown as solid elliptical curves. The region favoured by WMAP3 
data is located in the centre of the figure not overlapping with any of the the EDM allowed 
domains. 



to locate regions in which large phases are obtained satisfying all experimental bounds. The 
role of the phase ^3 is vital in this approach. For fixed .^2 one varies ^1 until the electron EDM 
de is small, by cancelling chargino against neutralino contributions. Subsequently .^3 is rotated 
until the neutron's EDM, dn, is made small. 

This is the prescription followed in [20] to make both de and dn lie within their experimental 
limits. However as discussed in the previous section due to the two-loop RGE dependence of 
de on ^3 the cancellation in dg may be lost and therefore ^1 need be re-rotated. It should be 
remarked that for large tan/3 the two-loop contributions are important and cannot be ignored. 
In those cases therefore we apply this recipe by taking into account the appearance of two-loop 
contributions, as well as additional contributions from other sources, such as chromoelectric 
dipole and gluonic dimension-6 operators, which may be important. Therefore following this 
procedure one may be able to find points in the Ci,2,3 parameter space yielding small de,dn. In 
this case however the Hg dipole moment is not guaranteed to be within its experimental limits. 
Besides, since the cosmologically allowed regions depend on these phases, even if we start from 
a point on which the relic density is acceptable and implement the cancellation mechanism, 
we do not necessarily end up with phases that respect the cosmological bounds put on the 
relic densities. Such a situation is depicted in Fig. [6l For the given inputs the phases Ci,2,3 sue 
rotated until the electron's and neutron's EDMs are within their experimental limits. This point 
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corresponds to .^2 = 0.007 vr and values of ^1^3 located at the centre of the small grey shaded 
region shown in the figure. This grey region is the overlap of the two stripes, one between 
the dashed lines, which designates the de acceptable region, and the other between the dotted- 
dashed lines which designates the region allowed by the neutron EDM bounds. Starting from 
this particular ^1,2,3 point, and keeping ,^2 fixed, we have plotted in the ^1^3 plane the allowed 
domains by EDM and relic densities. The allowed region by the Hg dipole moment is confined 
between the dotted lines (in magenta) which, as one can see, does not overlap with the grey 
region allowed by electron and neutron ED Ms. In this figure it is also seen how the phases are 
fine tuned, especially ^3, to achieve acceptable EDMs for the electron and neutron in the sense 
that only in a small region both df. and d„ can be within their experimental bounds. On the 
same figure the contours of constant neutralino relic density are shown as solid elliptical curves. 
In the example shown, one observes that the location of the region favoured by WMAP3 data 
occupies only a small portion at the centre of the figure at which ^1^2 — 0.0 not overlapping with 
the EDM allowed domains. This figure represents a rather typical example of the difficulty one 
encounters to reconcile both EDMs and cosmological bounds. 

In Fig. [7] we have taken moderate values mo = 500 GeV, M1/2 = 480 GeV, Aq = 100 GeV. 
and large tan /? = 50. The masses of the top and bottom are mt^b = 171.4, 4.25 GeV. The figure is 
constructed from one million random points in the ^1,2,3 parameter space. All other phases are as- 
sumed zero at the unification scale. Initially a particular point, ^1,2,3 = — 0.0727r, 0.9537r, 0.0477r 
is found by tuning the phases, in the way described earlier, satisfying all EDM bounds including 
in this case the EDM constraints from Hg atoms as well. The random sample is chosen to 
include this particular point. In Fig. [7] we classify these points in the planes ^1^3 (left panel) and 
^2,3 (right panel) according on what bounds each point satisfies. Light grey squares, forming the 
dense light grey region (region-1), includes points that satisfy only the Higgs and chargino mass 
bounds. The subset of these, shown as dark grey region (region-2), includes points that in addi- 
tion they satisfy the upper observational limit for the relic density, i.e. i^^^o < 0.117 but they 
do not necessarily fall within the region dictated by the WMAP3 region. Their subset, marked 
as triangles (region-3), are within the 2cr WMAP3 limits for CDM, 0.089 < Q^/i^ < 0.117 . For 
the rest ^x^q < 0.089 and therefore they are of relevance if additional components, except the 
LSP neutralino, contribute to the total DM density. 

The points shown as coloured-filled diamonds satisfy the electron and neutron EDM bounds. 
Of those only their subset filled with different colour (red) satisfy, in addition, the Hg EDM 
bounds. If any of these points falls on regions-1,2 or 3, defined before, then in addition it satisfies 
the corresponding bounds designating this particular region. One observes that irrespectively of 
the cosmological, and other accelerator data, the EDM bounds by themselves are hard to satisfy 
for moderate mcAf^rt. Despite the fact that a particular point was found which respects all 
three EDM bounds, the random sample of one million points in the ^1,2,3 space leaves only 
a few that observe the EDM limits for electron and neutron and even fewer that satisfy all 
three EDM constraints. This demonstrates that the values of the phases must be fine tuned 
to comply with experimental data of electric dipole moments. In the examples shown none of 
the diamond points satisfy the WMAP3 cosmological constraints. In fact for these points the 
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Figure 7: Scatter plots in ^3,^1 (left panel) and ^3,^2 (right panel) plane, based on a random 
sample of the gaugino phases ^1,2,3 for fixed mo = 500 GeV, M1/2 = 480 GeV, ^0 = 100 GeV 
and tan/3 = 50. The remaining phases at the unification scale are taken zero. The light grey 
squares represent points that satisfy the mass bounds for the light Higgs boson and other SUSY 
particles. The dark grey region is formed from points that in addition they satisfy the upper 
WMAP3 bound for the relic density f^-^/iQ < 0.117. The yellow triangles represent their subset 
falling within the WMAP3 region, 0.089 < ^x^q < 0.117. The one-colour (green) diamonds 
satisfy the electron and neutron EDM bounds, while those filled with different colour (red), 
satisfy, in addition, the Hg EDM bounds. Depending on which region they lie on, they may 
observe the bounds put by WMAP3. 
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Figure 8: Scatter plots in ^i,Mi/2 (left panel) and ^3,Mi/2 (right panel) plane, for a random 
sample with ^2 fixed and random values for ^1^3 and M^^. The notation is as in Fig. [71 

predicted neutralino relic density is below the limits put by WMAP3 and therefore additional 
DM candidates must exist to fill the deficit. 
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In Fig. [8l with the same starting values for the phases ^1,2,3 and the parameters ttiq, -^1/2 > 
Aq, tan/3, for which agreement with all EDM bounds are obeyed, we keep ^2 fixed and generate 
a random sample of one million points for ^1^3 and M^ rt. The sample includes the starting values 
for ^i^3,Mi/2- In Ci!-^i/2 (Isft panel) and C3,Mi^2 (right panel) planes, the points satisfying the 
particular criteria as described in Fig. [7| are displayed. The notation is as in Fig. [71 The points 
satisfying the electron and neutron EDM bounds are few, and only one satisfies all three EDM 
bounds. In the case considered a few of these points overlap with the region-3 (yellow triangles) 
and therefore for these all available data are obeyed, with the exception of the Hg EDM bound. 
The limited number of points satisfying the EDM bounds, in this case too, it indicates that the 
phases must be fine tuned to agree with experimental data. 

The Hg EDM bound, in general, poses a severe obstacle in obtaining agreement with cos- 
mological data. However by tuning appropriately the gaugino phases there are cases where 
agreement with all ED Ms is obtained at a particular point tjiq, Mii2-,Aq of the parameter space. 
Then by varying mo,Mi/2 there is a chance that one succeeds in obtaining cosmologically and 
EDM allowed regions which overlap. The chance of obtaining this is increased if one considers 
funnel regions of rapid neutralino annihilation via a Higgs resonance, which occupy extended 
regions in the parameter space. In figures [9l and 1101 we display cases for low and large tan/3. The 
values for the phases ^1,2,3, in each figure, are fine tuned for specific mQ,Mi/2, around 600 GeV, 
in order to obtain agreement with all three EDMs, de, dn, dng, for the electron, neutron and 
Hg respectively. We have taken Aq = 100 GeV and the phase of /i is taken vanishing. All 
other inputs are shown on the figures. The dashed lines (in blue) delineate the boundaries of 
de, the dotted lines (in magenta) those of dng and the dotted-dashed lines (in orange) those of 
dn- If only one boundary line is shown it simply means the other one lies outside the displayed 
mo,Mi/2 range. The allowed by EDM constraints regions are between the boundaries in each 
case. In all cases displayed there are regions where all electric dipole bounds are simultaneously 
satisfied. Also shown are the Higgs line 114.0 GeV, long-dashed (red), and the line 115.0 GeV, 
long dashed-dotted (in blue), lying on the left almost vertical to the M1/2 axis. Their upper 
ends touch the no-electroweak symmetry breaking region, designated as a hatched area. At 
the bottom the triangle-shaped shaded region is excluded since there the stau is lighter than 
any of the neutralinos. The 2 — a cosmologically allowed WMAP3 regions are shown as shaded 
contours (in light green). 

On the left panel of Fig. [9] the value of tan /3 = 10 is small and the cosmologically allowed 
regions are not so extended. Although for the specific values of the gaugino phases there are 
large portions of the parameter space compatible with all EDMs, the overlap of these regions 
with the regions allowed by the WMAP3 data extend to large values of mo > 5.0 TeV, not 
shown in the figure, along the focus point region [84]. On the right panel of the same figure the 
value of tan /3 = 30 is considerably larger and the cosmologically allowed domains have a larger 
overlap with the regions allowed by EDMs. In this case there is a region in which all data are 
satisfied for tuq ~ 1 TeV and M1/2 > 2.2 TeV due to the fact that the allowed by dg domain 
is enlarged, allowing for smaller M1/2 values, which includes part of a funnel that starts being 
formed. This is located on the right of the figure, just above the shaded region which is excluded 
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Figure 9: The mo,Mi/2 parameter space for non zero gaugino phases, with values shown on 
the figures, ^o = 100 GeV and the phase of fi is taken vanishing. The dashed hnes (in blue) 
delineate the boundaries of de , the dotted lines (in magenta) those of dng and the dotted-dashed 
lines (in orange) those of d„. If only one boundary line is shown it means the other one lies 
outside the displayed mo, M1/2 range. The allowed by EDM constraints regions are between the 
boundaries in each case. The rest of curves and shaded regions are as in Fig. [5j 
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Figure 10: As in Fig. [9] for different values of tan/3 and the gaugino phases. 

since there the stau is the LSP. Due to the heaviness of M1/2 this is outside the reach of LHC. 
As in the left panel case there is also a focus point region, acceptable by all data, starting now 
from smaller values of rriQ > 2.1 TeV, which tracks the border of the no-electroweak breaking 
region. Part of it includes points with M^ « < 800 GeV being therefore accessible to LHC. 

Notice that the values of the phases for the right panel are different from those of the case 
displayed on the left. Due to the fine tuning of the phases chosen increasing the value of tan/? 
from 10 to 30 in the left panel we do not get a picture resembling the one we display on the 
right panel. 
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Figure 11: As in Fig. [H with a non- vanishing phase of /i and values of the gaugino phases and 
tan P shown in the panels. 

In Fig. [lOl for different sets of the gaugino phases and large values of tan/3, we display 
regions that are allowed by all data. Again the phases have been fine-tuned and are different 
for each case shown in the two panels. On the left panel the value of tan/3 is large, tan/3 = 45, 
and the cosmologically allowed domains are funnel-shaped, extended diagonally towards high 
77^0 and Ml/2 values, tuq ~ 3.0 TeV, M1/2 — 2.2 TeV. At the same time the acceptable EDM 
domains are also extended having a large overlap with the cosmologically allowed region. In 
this particular case there is a large portion of the parameter space for mo,Mx/2 > 1-3 TeV 
in which all experimental data are satisfied. A small part of this region is within the reach 
of LHC. The conclusion from this is that by tuning the gaugino phases there can be found 
extended regions in the parameter space in which rapid neutralino annihilation through a Higgs 
resonance can coexist with regions allowed by the stringent constraints imposed by the electric 
dipole moments. On the right panel the value of tan/3 = 50 is larger and it also allows for the 
formation of rather extended cosmologically allowed regions which are almost funnel-shaped, but 
no so peaked as the case considered previously. The boundaries of EDMs in the case shown lie 
within the range mo < 2.0 TeV, M1/2 < 1.0 TeV and one observes that EDM and cosmologically 
allowed domains again overlap. The overlapping regions are not so extended however, as in the 
case considered previously. In the case shown it is only a small region centered around the point 
mo,Mi/2 = 900,600 GeV. One observes that larger values of tan/?, unlike the CP-conserving 
case, do not always correspond to funnels spanning the highest m,o,i/2 regions. This is due to 
the sensitivity of the bottom mass, and hence the funnel regions, with the CP-violating phases. 
As an effect cosmologically allowed funnel-shaped regions can appear for smaller values of tan /?, 
as compared with the CP-conserving case, and can coexist with EDM allowed domains. 

In the previous analysis the n phase has been put to zero but one can also seek for cases 
where large values of this phase are allowed. This is forbidden in cases where the gaugino phases 
are switched off. In mSUGRA models, in which the common phase of the trilinear coupling 
and the /i phase are the only allowed super symmetric CP-violating sources, it is known that the 
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Figure 12: Scatter plots based on a random sample with mo = 2000 GeV, M1/2 = 1720 GeV, 
^0 = 100 GeV, tan/3 = 45 and the phases in the region -0.87 < ^i/vr < -0.77, 0.34 < ^2/vr < 
0.44, —0.46 < ^s/vr < —0.36 and 0.55 < (P^/tt < 0.65. The dots (in red) are points of the random 
sample that satisfy all the EDM and cosmological constraints. The (blue) diamond marks the 
values of the phases used in the right panel of Fig. [TTl 

phase of fi is tightly constrained by the EDM data, especially by this of the electron. However 
when one allows for the presence of different gaugino phases the situation is altered and this 
restriction is relaxed. In Fig. [TT] we display a case where (f)^ is non- vanishing and large. On the 
left panel a case with tan/? = 10 is shown. On the right panel tan/? = 45. The value of Aq is 
Aq = 100 GeV. In both cases the phase of fi is non-vanishing and there are regions compatible 
with EDMs, cosmological data, as well as all other accelerator data. On the left panel only the 
focus point region is compatible with all data, while on the right panel both focus point region 
and the cosmologically allowed funnels are allowed. Both regions include points accessible to 
LHC, which are larger as compared to the cases considered previously. 

Note that the values of arg{fiMi^2) at Mqut are sizable, 0(0.1 vr). Due to absence of 
renormalization of the ^ phase with the energy scale and the small renormalization of the gaugino 
phases, which at one-loop are independent of the SUSY inputs at Mqut, these quantities retain 
almost the same values at the EW scale Mz, for every point in the m0)-^i/2 plane. The 
difference of their values at Mz from those at Mqut is at per mille level which is very small and 
important only for EDMs. As stated in the introduction these combinations of phases set the 
measure of sufficient Higgsino and gaugino driven Baryogenesis at the EW phase transition. On 
these grounds and in conjunction with the fact that the magnitude of ^ is comparable to that 
of Mi^2 in a large region of the allowed parameter space, displayed in this figure, these regions 
may be also compatible with these Baryogenesis mechanisms [49-51]. 

In order to show how sensitive the selected regions are to variations of the phases, which 
were chosen to suppress the EDMs via the cancellation mechanism, we pick a particular point 
located in the allowed funnel region on the left panel of Fig. [TTJ Then by varying the phases 
Ci,2,3 and (p^ around the values displayed in this figure, we produce a random sample consisting 
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of 100,000 points. The scatter plots displayed in Fig. [T2l are based on this random sample 
with mo = 2000 GeV, M1/2 = 1720 GeV, ^0 = 100 GeV, tan/? = 45 and phases in the region 
-0.87 < 6/7r < -0.77, 0.34 < ^/^r < 0.44, -0.46 < 6A < -0.36 and 0.55 < (/>/,/7r < 0.65. 
The scattered dots ( in red ) represent the random points that satisfy all the EDM bounds and 
the WMAP3 bound on the neutralino relic density. The (blue) diamond in the center of each 
panel, marks the point corresponding to the values of phases on the right panel of Fig. [TTJ One 
can see that the allowed variations on 0^ and ^3 are of the order of ~ 0.057r rad. The same applies 
to the phase ^2 as well. On the other hand, the range of the values of the phase ^1, that are 
compatible with the EDM and cosmological constraints, appears to be much broader. Therefore, 
in the large tan /3 regime, the amount of tuning required to locate extended cosmological funnels 
that are also compatible with the EDM bounds, is of the order of O.Olvr rad for ^2,3 and t/)^ and 
O.lvr rad for ^1. In the displayed results the phases of the trilinear couplings at the GUT scale 
are taken zero. However we have checked numerically that variations on them of the order O.lvr 
rad do not destabilize the cancellation mechanism for the EDMs and therefore the fine tuning 
of these phases is less restrictive. 

6 Conclusions 

We have considered supersymmetric models in the presence of CP-violating sources residing in 
the Higgsino-mixing mass term and SUSY breaking soft parameters. In the simple, mSUGRA 
based, versions of such models, that have been extensively studied in the past, there are two 
independent phases at the unification scale usually chosen to be the phase of the Higgsino mixing 
parameter /x, 0^, and the phase of the common trilinear coupling ^O) 4' A- The application of the 
EDM and cosmological bounds constrains the phase of /i to be (/>^/vr < 0.01 rad, whose small- 
ness poses severe obstacles in certain Baryogenesis mechanisms, while (J)a remains practically 
unconstrained and can be at least ten times larger. 

In our study we have used the revised top mass, whose central value is mt = 171.4 GeV, and 
scan the parameter space allowing for non-universal boundary conditions for the phases at the 
GUT scale. Such models are described by thirteen CP-violating phases, residing in the gaugino 
mass terms, the trilinear scalar couplings and the /i parameter, one of which can be rotated away. 
An additional phase which misaligns the VEVs of the Higgs doublets is determined from the 
minimization conditions and affects the analysis. We follow a top-bottom approach according to 
which the low energy values of all parameters involved, including their magnitudes and phases, 
are determined from their values at the unification scale after running of the renormalization 
group equations. In such an analysis we showed that the two-loop running of the CP-violating 
phases induces important corrections to the electric dipole moments of the known species, which 
are absent in the one-loop running analysis. This results to further constraints of the allowed 
parameter space. Important corrections to EDMs are also induced by the one-loop misalignment 
angle of the Higgs vacuum expectation values which are augmented for large tan (i and in the 
presence of a non-vanishing gluino phase. 

The role of the gaugino phases is particularly important allowing for suppression of electric 
dipole moments if they are fine tuned at the unification scale. By tuning appropriately the 
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gaugino phases ^^1,2,3 and the phase (j)^ of /i at the unification scale, there can be found regions 
in the parameter space, along the focus point and neutralino pair annihilation through a Higgs 
resonance, which are permitted by dark matter WMAP3 cosmological constraints, and are also 
compatible with the electric dipole moment limits and all accelerator data. 

This can be accomplished for large values for the phases of the gaugino masses and the ^ 
parameter, of the order of the 0(1) vr rad, which are however fine-tuned. The phases ^2,3 and 
(j)^ need to be adjusted at the O.Olvr rad level, while the adjustment of ^1 is by an order of 
magnitude less restrictive. The phases of the trilinear scalar couplings although they play a 
key role for EDMs are not fine tuned when the magnitude of the trilinear coupling is small. 
The role of the gluino phase is important in this analysis. It can be used, by implementing 
the cancellation mechanism, to make neutron's electric dipole moment small, and in addition it 
allows for a large non-vanishing phase for /i. Large values of ^ and gaugino phases are needed 
in certain electroweak Baryogenesis mechanisms . Higgsino and gaugino driven Baryogenesis 
requires the magnitude of /i to be comparable to the gaugino masses which can be fulfilled in 
the major portion of the parameter space in the constrained models considered in this work. 

Since the suppression of the EDM is feasible for small and intermediate values of the SUSY 
breaking scale, there are regions in the parameter space which are within the reach of the 
forthcoming LHC experiments. 
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1 Appendix 

In this section we present the basic expressions and conventions used in this work so that a direct 
comparison with other works and different notations is made possible. The supersymmetric part 
of the Lagrangian is specified by a superpotential given by 

W = htQ^ eH2U^ + h Q^ eHiD^ + K L^ e H2 E^ + fi Hf e H2 , (A.l) 

where the elements of the antisymmetric 2x2 matrix e are given by ei2 = —£21 = 1 . In the 
superpotential above we have only shown the dominant Yukawa terms of the third generation 
and we do not allow for flavour mixings. According to this the Yukawa couplings of the left. 
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right handed top and bottom niultiplets are 

W = htH^UU^-hb Hi D D^ 
The scalar soft part of the Lagrangian is given by 

'-'Scalar ^ / ^ ^^i I 'Pi 



|2 



- {At ht Q^ eHiU^ + Ab h Q'^ e Hi D'' + A^ K L^ e H2 E" + h.c.) 

- ( ml Hf e H2 + h.c.) , (A.2) 

where the index i in the sum in the equation above runs over aU scalar fields and all fields 
appearing denote scalar parts of the supermultiplets involved. 
The gaugino fields soft mass terms are given by 

Cgaugino = -^ { Mi B B + M2 w'-'^ VF^^) + M3 G G + h.c.) (A.3) 

In this equation B, W^^\ G are the gauge fermions corresponding to the U{1), SU{2), 
and SU{3) gauge groups. Our notation is that of Ellis and Zwirner [86] with the signs of the 
gaugino masses and that of the bottom and tau Yukawa couplings reverted. The one- loop RGEs 
found in that reference coincide with ours if the signs of the gaugino masses are flipped. RGEs 
of course are insensitive to the Yukawa sign convention. 

The two-loop renormalization group equations (RGEs) in the general case, including super- 
symmetric CP violations, can be read from [72] and [73]. In [72] the soft SUSY breaking part 
of the Lagrangian is written as 

£ = _1 h'^k ^^^^^^ _ ^f^ij ^.^. _ 1 ^rnyiM*^ - ^M(^a)XaXa + h.c. (A.4) 

Keeping only flavour diagonal mass terms for the scalars (pi, as we have assumed throughout, 
the mass terms become 



/I 



Prom Eq. ([T]) we observe that the imaginary parts of {rn )^ do not appear in the Lagrangian. 
Besides it is easy to see that the RGEs of the quantities appearing in the Lagrangian do not 
depend on the imaginary parts of {'m'^)^. Therefore they do not affect the physical quantities 
and can be taken anything. 

Por the two-loop RGEs, those of the trilinear couplings At^h,T are found if one replaces hijk 
and the Yukawas couplings Yijk in the RGEs found in [72] by /i = —A ■ Y, where we have 
suppressed the flavour indices, and identifying the gaugino masses used in that reference with 
ours. The one-loop parts of the RGEs we get with such a replacement coincide with our one- 
loop RGEs for the A^^i^t which are identical, modulo the sign difference in the gaugino masses 
as mentioned before, with reference [86]. Note that the same would not have happened have 
we used h = —A ■ Y. To complete the correspondence and having as our guideline the one-loop 
results for the RGEs, fj, of [72] should be replaced by our fi and the coupling B of the two 
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Higgs scalars should be replaced by our —m^. With this correspondence the two- loop RGEs are 
retrieved unambiguously from [72] and are adapted to our notation. 

Concerning the neutralino and chargino mass matrices, In the B, W^'^' , iH^, iH2j basis the 
neutralino mass matrix is 



M 



N 



/ Ml g'{Ht) -g'{m)\ 

M2 -9 {HI) g{H*2) 

g'{H*) -g{H*) -/. 

\-g'{m) g{H*2) -f, / 



(A.5) 



The field dependent neutralino mass matrix entering into the effective potential has a similar 
form with the Higgs VEVs replaced by the neutral components of the corresponding Higgs 
doublets. Note that the gaugino masses Mi^2 the parameter^ as well as the Higgss VEVs are 
complex in general. Regarding the Higgses one can verify that the mass eigenstates depend only 
on their relative phase. One can diagonalize the symmetric neutralino mass matrix as 

O'^MnO =diag(m^o,m^o,m^o,m^o] , (A.6) 



where O is complex in general. In our approach this matrix is chosen is such a way that the 

eigenvalues m^o are real and positive. 

The chargino mass matrix can be obtained from the Lagrangian terms 



^ThZUos = - {W^,^H^) M, ( ,^+ ) + (h.c) , (A.7) 

where we have defined W = — — ^ii^- — . The mass matrix is 

v2 

( M2 - g (m) \ 
Mc = { ' , A.8 

V - ^ ^^1) 1^ J 

which can be diagonalized as 

UMcV^ = ("^^^ M . (A.9) 

Thus, 

^ThVrginos = -m^.XlXl - m^2X2X2 ■ (A.IO) 

The Dirac chargino states xi,2 in this equation are defined by 

with the two component Weyl spinors A^^ 2 related to W , iH^ , iH2 by 



V 



m 



) ^ ( 4) ' (^ ' '^' ) ^^ ^ ^^' ' ^' ^ ■ ^^'^^^ 
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The gauge interactions of charginos and neutralinos can be read from the Lagrangiaqfl 

-C = 5 {W;r + W^Jt) + eA^Jl^^ + -Z^J^ . (A.13) 



(A.14) 



e 
sc 

In the equation above s = sin^v^! c = cos 9w ■ Also, 

The currents J^, Jem and J^ are given by 

Jt^Xar[rLr!:, + VRV^]Xi a = l..A, i = l,2, (A.15) 

where Vl^r = -^^ and 

r^^+^Ol^V^^-O^^V*, , v!^,^-^0saU:,-02aU:,. (A.16) 

The electromagnetic current Jl^m is 

JL = XiT'^Xi + X27^X2 . (A.17) 

Finally, the neutral current J^ can be read from 

4 = hr [ VlAJ- + VnAf^ ]x, + \ Xal^ [ 'PLB^b + -PrB^, ] xl , (A.18) 

with 



A^- = C%--F.2^V 



Aj = C 6ij - -Ui2U*2, 



2 



,L 1 



B^b = -B':C- (A.19) 

The chargino and neutralino couplings to sfermions are given by the following Lagrangian 
terms 

C- = iXi {Vl a{/ + Vr b{/) f /; + i ^, {Vl 4[ + 'Pr H[ ) f fj + (h.c) . (A.20) 

In this, Xi {i = 1;2) are the positively charged charginos and xf the corresponding charge 
conjugate states having opposite charge. / , /' denote "up" and "down" fermions, quarks or 
leptons, while /j , f[ are the corresponding sfermion mass eigenstates. The left and right-handed 
couplings appearing above are given by 

a{f = gU,iKf[ + hf,Ui2Kf^ , Hf = hfVi2Kf[. 



In our notation e =electron's charge. 
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In the equation above hf , hf/ are the Yukawa couphngs of the up and down fermions respectively. 
The matrices K^'-l' which diagonaUze the sfermion mass matrices become the unit matrices in 
the absence of left-right sfermion mixings. For the electron and muon family the lepton masses 
are taken to be vanishing in the case that mixings do not occur. In addition the right-handed 
couplings, are zero. 
The corresponding neutralino couplings are given by 

C = i Xai-PL «{/ + Vr 6{/) ff* +iXa {Vl a{/' + Vr b{/') f !'* + {h.c) . (A.21) 

The left and right-handed couplings for the up fermions, sfermions are given by 

ail = V2{gTp2a+g'^Ola)Kf^ + /l/ 04a i^/2 , 



while those for the down fermions and sfermions are given by 

JJ = V2{gT],02a+9'^Oia)Kf[ - hf.O^aKf^ 
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